Effective models of two-flavor QCD: flnite and mg-dependence 



PACS numbers: 12.38 Aw, 12.38 Mh 



T. Kahard 

Department of Physics, University of Jyvaskyla, P.O.Box 35, FIN-40014 Jyvdskyla, Finland 
and Helsinki Institute of Physics, P.O.Box 64, FIN-00014 University of Helsinki, Finland 

K. Tuominen^ 

CP^ -Origins, Campusvej 55, DK-5230 Odense M, Denmark 
and Helsinki Institute of Physics, P.O.Box 64, FIN-00014 University of Helsinki, Finland 

We study effective models of chiral fields and Polyakov loop expected to describe the dynamics 
responsible for the phase structure of two-flavor QCD at finite temperature and density. We consider 
chiral sector described either using linear sigma model or Nambu-Jona-Lasinio model and study the 
I phase diagram and determine the location of the critical point as a function of the explicit chiral 

symmetry breaking (i.e. the bare quark mass rrig). We also discuss the possible emergence of the 
quarkyonic phase in this model. 
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I. INTRODUCTION 

O I Hadronic matter undergoes a phase transition from hadronic matter into a partonic matter at high temperatures or 
JL densities. To predict the equation of state, study the existence of possible critical point(s) in the (T, /i)- phase diagram 
and the properties of the phase transitions presents a theoretical challenge on studies based on the fundamental 
theory of strong interactions, QCD. To obtain some insight into the QCD dynamics of quarks in the nonperturbative 
'— ' domains, models like the Nambu-Jona-Lasinio (NJL) model have been developed. These models are based on the 

^ chiral symmetry of light quarks. 
^ Based on generic effective theory methods, identification of relevant degrees of freedom and symmetries, effective 

^-H models able to account for the two most important features of QCD, deconfinement and chiral symmetry breaking, 
fT^ have been developed and studied [3 [5] . These models include both chiral fields and the Polyakov loop as relevant 
On degrees of freedom. The chiral symmetry is effectively represented by a NJL model or linear sigma model (LSM), 
while the Z3 symmetry relevant for deconfinement in pure gauge theory is described by a mean field potential for 
\^ Polyakov loop. The key role is played by the interactions coupling these two sectors. In the limit of light quarks the 
dynamics are driven by the chiral degrees of freedom and the decrease in the chiral condensate as the temperature is 
increased results in increase of the Polyakov loop which in turn results in deconfinement and explains the intertwining 
of these two seemingly unrelated features to a single phase transition fT-'T] . On the other hand, if the current quark 
^ masses are taken large, chiral symmetry broken explicitly, the dynamics is close to that of pure gauge theory, i.e. 
k> dominated by the Polyakov loop. Nevertheless, at deconfimenent the interactions now lead to decrease in the chiral 
rN condensate and the two transitions again coincide. In real QCD approximate chiral symmetry is typically expected 
to play the dominant role. 

In an earlier work we have compared the models of two-flavor QCD where the chiral sector is represented either 
with LSM or NJL model [5] at finite temperature and density. We found that while at zero chemical potential both 
cases give practically coincident results, the relative uncertainties increase when finite quark densities are considered. 
In particular the predictions for the location of the possible critical point in the (T, /x)-plane differ widely [8]. As 
a further direction to study and constrain these models we considered the dependence on explicit chiral symmetry 
breaking, i.e. the value of the quark mass m,, in these models [5]. In this brief report we complete this previous study 
by extending the analysis to finite density. We briefly comment also on the possible emergence of a quarkyonic phase 
introduced in recent literature jlOj. 

The paper is organized as follows: in section ITT] we briefly recall the basic definitions of the two models we consider 



and explain the model parameters which allow for arbitrary pion mass. In section III we present our main results and 
in section Hvl our conclusions and outlook. 



t On leave of absence from Department of Physics, University of Jyvaskyla 

* Electronic address; topi.kahara@ j yu^ 

-^Electronic address: kimmo.tuominen@phys.jyu.fi 



2 



II. THE MODELS 



A. Models at the physical point « 140 MeV 

As explained in the introduction, in this work we continue our study 9J of the quark (or pion) mass dependence 
of the QCD phase diagram. We consider the PNJL and PLSM models which consist of a chiral part, a lattice fitted 
Polyakov potential and a simple interaction between the two. The study of these models is done in the mean field 
approximation. For a detailed description of the derivation see [9], here we will simply state the resulting grand 
potential 



^ = t^chiral + Uf + i}q 



(1) 



The two models differ only in the chiral part which corresponds either to the NJL or LSM models. The linear sigma 
model (LSM) consist of the sigma meson and the pions with their mutual interactions, and interactions with quarks. 
The NJL model on the other hand describes only quarks with an effective four-fermion interaction. In both models, 
the explicit chiral symmetry breaking is taken into account. 

The deconfining phase transition is included in both models through the mean field potential 



(2) 



where 



62 (T) ^ ao + ai 



a2 



(3) 



and the constants 0^,6^ are fixed to reproduce pure gauge theory thermodynamics with phase transition at Tq = 270 
MeV; We adopt the values determined in |4j , and shown for completeness in table |l] Here i is the gauge invariant 
Polyakov loop in the fundamental representation. Instead of the polynomial form for the Polyakov loop potential 
([2]), also other possibilities exist. An example is the one introduced in [11] which has the advantage that £ is always 
confined to the values between zero and one. However as discussed in [1^ these two possible forms for the potential 
do not differ significantly at temperatures below 300 MeV, the region we are interested in. Further improvement to 
the potential would be the inclusion of /^-dependence in the Tq parameter as is done for example in fl3j and [14j : we 
will discuss this briefly in Sec. |I1I C 



The chiral potentials in ([T]) are 



^chiral 









(7) 
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HM 



for LSM 



chiral 



2G 



for NJL. 



(4) 



(5) 



At the physical pion mass the parameters in the above equations are fixed by the physical vacuum properties. In the 



LSM model H = f^ml and = fl 
determined by the tree level mass 



m^/A^, where = 93 MeV and ra^ = 138 MeV. The coupling ss 20 is 
= 2A^/^ + ml, which is set to be 600 MeV. In the NJL model we fix the bare 
quark mass to be ruq = 5.5 MeV and the coupling G = 10.08 GeV~^. The constituent masses M are related to the 
qq and a expectation values throught the relations M — niq — G{qq) in NJL and M = g{a) in LSM. In the latter case 
the coupling constant g is fixed to 3.3 corresponding to the baryon mass ~ 1 GeV. 

The final term in ([T]) includes the interaction between the chiral and Polyakov sectors and reads 



^qq 



-2NfT 



(27r)3 



Tr^ In 



l + Le-(^-^)/^ 



Tr^ In 



(6) 



'gsAa/T]. The above contribution is of the same basic form for both 



where the Polyakov loop matrix is L — exp[ 
PLSM and PNJL models with E = i/p 2 + ][.p, where M is the constituent mass of the model in question as defined 
above. The trace over color remains and using the definition of the Polyakov loop £ — (Trc(£))/iVc and taking the 
Polyakov loop matrix L corresponding to a static background field Aq, one obtains 



= -2NfT 



(2^ 



(in [1 + 3(£ + r e-(^-'^'/^)e-(^-^)/^ + e-3(^-^)/^' 



+ In [1 + 3(r + le-^E+f.)/T-^^~(E+t.)/T ^ g-3(B+M)/Tj j 



(7) 
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The interaction potential flqq includes also a vacuum term omitted from equations ^ and Q 

~6Nj I ^Ee{A'^\p\'). (8) 

This term is neglected in the PLSM model, but included in the PNJL model where it is controlled by the cut-ofF A, 
which we set at 651 MeV. A summary of the parameters at the physical point is shown in table |l] 

The thermodynamics of the models are determined by solving the equations of motion for the order parameters, 

8M='^ aZ=«' 

and then the pressure is given by evaluating the potential on the minimum, p = —il{T,fi). We have chosen the 
constituent mass M as a basic variable since this most conveniently allows us to discuss both models simultaneously. 
It is also straightforward to write the results in terms of the condensates (cr) and (qq) since these are linearly related 
to M in each case. At finite chemical potential the mean field potential SI is complex due to the Polyakov loops and 
minimizing such a potential is meaningless. A simple way to overcome this problem is to treat the Polyakov loop 
parameters £ and £* as independent real variables, which will be a sufficient approximation for our current analysis. 
Inaccuracies of this treatment as well as improved methods have been discussed for example in O lllj . 



TABLE I: The parameters used for the effective potential 



LSM: m„ 

93 MeV 138 MeV 600 MeV 

g \ H 

3^3 4.44 1.77- 10"^ GeV^ 

JUL: nTq A G 

5.5 MeV 651 MeV 10.08 (GeV)"^ 

Polyakov: ao ai 02 

6.75 -1.95 2.625 

as &3 64 

-7.44 0.75 7.5 



B. Models away from the physical point 

To explore the region of parameter space in which the pion mass differs from its physical value, we need a consistent 
way of setting the model parameters in a such region. In the PNJL model this is quite easy since the bare quark mass 
is a direct input parameter that controls the amount of chiral symmetry breaking in the model and the pion mass 
is then calculated from the model. In the PLSM case the situation is not so simple since we have four parameters, 
m^, rria- and g, which arc connected with each other. To make the comparison with the PNJL model easier we 
introduced in [5] a parametrization based on lattice results P"5HT7] which relates the above parameters directly to the 
bare quark mass ruq. Here we will give a brief summary of this parametrization. 

The lattice data connects the pion decay constant, pion mass and sigma mass through the equations 

mla'^ = {Ai{mqa)^ + B{mqaf) (10) 
V2Ua = 0.06672 + 0.221820 x (mqa) - V2Ca (11) 
m, = ^ml + D (12) 

obtained from ,15J and [16^. The parameters in the above equations are shown in table [lT| and chosen to reproduce 
the physical vacuum values for the PLSM model parameters i.e. f-^ = 93 MeV, = 138 MeV and rria = 600 MeV 
when the quark mass rriq is set to 5 MeV. 

The PLSM model parameter g is determined through the relation 17 /jr = Mn/S, where the nucleon mass Mjf is 
parametrized in the form 

Mm = Mo + 4Cim2 - -|%m3 . (13) 
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TABLE 11: The lattice parameters 



a Ai S B C 

0.505306 (GeV)"^ 0.82725 0.16413 1.88687 1.18 MeV 

D ^ Mo Cl QA 

565.15 MeV 1.83 (GeV)-^ 868 MeV 0.9 (GeV)-^ 1.267 



This is a chiral perturbation theory fit truncated to 0{ra^). In (T7j it has been shown that for a good description of 
nucleon mass one should keep terms up to and including ©(m^), but for simplicity we have chosen the truncated fit. 
Previously, in [5], we used the same nucleon mass formula truncated to ©(m^) but this exaggerated the strength of the 
coupling g at larger vriq to an extent that was found to have a large effect on the phase diagram; we chose to improve 
by adding the C'(m^)-term. The parameters are chosen so that the PLSM model vacuum values are reproduced and 
the corresponding parameter values are shown in table |lTj 

A comparison between the PNJL model and the PLSM model, now equipped with our lattice based parameter fit, 
is shown in Figure [l] Both the pion masses and the pion decay constants agree very well between the two models 
at low quark masses. As the bare quark mass, i.e. the amount of explicit chiral symmetry breaking is increased, the 
models start to deviate. At niq = 250 MeV, the largest quark mass shown in Figure fll the deviation is around 20 % 
for both the pion masses and decay constants and keeps increasing as one increases oare quark mass further. This 
growing deviation is a natural indication that the models, based on approximate chiral symmetry, start to fail as the 
explicit chiral symmetry breaking becomes large. It should be noted that the comparison between the PNJL and 
PLSM models includes no tuning of the coupling G or the cut-off A of the PNJL model. Actually we have checked 
numerically that the agreement between the PNJL and the lattice fitted PLSM pion mass curves in Figure [l] cannot 
be improved by altering the values of G or A from those shown in Table |l) 

In the following section we will study mainly the effect of explicit chiral symmetry breaking on the thermodynamics 
of the models. We focus in particular on the (T, /i)- phase diagram including possible critical points. 




FIG. 1: Pion mass and pion decay constants /.^ as functions of the bare quark mass in the PLSM and PNJL models. 



III. RESULTS 



By construction, in these models there are a priori two transitions: The chiral transition due to the (approximate) 
restoration of chiral symmetry and the deconfinement transition encoded into the Polyakov potential. The transitions 
can be studied through their respective order parameters, the constituent quark mass M and the thermal average 
of the Polyakov field, i. To both transitions one can assign their own critical temperatures. The definition of the 
critical temperature, however, is vague especially in the regions where the transition is a crossover and the order 
parameter shifts continuously. Since this is the case over large portion of the (T, /i)-plane, in this work we primarily 
define the transition temperature as the temperature at which the temperature derivative of the order parameter has 
a maximum. Even this definition has some problems, since in some cases the derivative has several local maxima 
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TABLE III: Bare quark masses with corresponding pion masses and critical points for both models. All values are in MeV. 



PLSM 






PNJL 






niq 




CP (T,m) 


ruq 




CP (T,^) 


0.1 


26 


None 


0.1 


19 


(147, 270) 


2 


93 


None 


2 


85 


(111, 306) 


5 


138 


(196, 135) 


5.5 


140 


( 88, 329) 


15 


222 


(168, 235) 


15 


231 


( 59, 364) 


50 


377 


(120, 353) 


50 


421 


( 58, 435) 


100 


518 


( 30, 455) 


100 


603 


( 92, 496) 


150 


629 


None 


150 


752 


(120, 536) 



indicating rapid changes at several different temperatures. The critical temperature is identified with the maximum 
at which the change in the absolute value of the corresponding order parameter is largest. Alternatively one could 
use the susceptibilities to define the critical temperature. 



A. The chiral phase diagram and the critical point 



The chiral transition can be determined by finding the temperature corresponding to the fastest change in the 
constituent quark mass M at fixed chemical potential (or vice versa). This transition temperature corresponds in 
most cases to the temperature at which the constituent mass drops below 50% of its vacuum value, only at large 
and fjL does the fastest change occur at a different temperature than the one where the decrease in the absolute value 
of the constituent mass takes place. Figure [2] shows the chiral transition lines in the (T, /i)-plane for different quark 
masses for both models; also the critical points are shown. The critical points indicate the points where a line of first 
order (discontinuous) transitions ends and turns into a crossover (continuous). 

As seen in Figure [2] the qualitative features of the phase diagrams in the two models are very similar: As the quark 
mass rises, the area under the transition line expands and the critical point moves towards larger /z. The quantitative 
difference in the transition temperature between the models is below 15% for the shown quark masses. However the 
critical points appear to be located quite differently in the models as we already noted for the physical value of 
in [5] . In the PLSM the transition for the lowest quark masses shown in Figure [2] is first order all the way so there is 
no critical point in the (T, /i)-plane, the same holds for the largest quark mass but now the transition is a crossover 
throughout the plane. Therefore, in the PLSM model the critical point is present only at a finite nuq interval and 
outside this interval the transition is either entirely crossover or entirely of first order. In |18j a qualitatively similar 
result has been obtained for a three flavor linear sigma model. The fact that the transition in PLSM is of first order 
over the entire (T, /x)-plane in chiral limit is due to neglect of the fermion vacuum energy |19| . 

In the PNJL case, where the fermion vacuum contribution is included, the critical points persist even for the smallest 
quark masses shown and the transitions at zero chemical potential remain crossovers. Also the critical temperature 
at the critical point starts to rise again at larger quark masses and does not disappear as in the PLSM case. A similar 
effect has been observed in pT, where a saturation of the critical point temperature was mentioned and attributed 
to a diquark dominated phase. Since our work does not include diquark degrees of freedom, we conclude that the 
behaviour of the critical point at large quark masses is a more generic feature of the PNJL model. 

It has been suggested in [20] that there might be, especially at small quark masses, multiple critical points in the 
(T, yu)-plane. We, however, found no evidence in either model to suggest that this is the case. In our previous work 
[9] we noted that at large quark masses the PLSM transition was first order at = giving some credence to the idea 
of multiple critical points at large quark masses. As mentioned in the previous paragraph and evident from Figure 
[2j this is not the case in our present work and the first order transition at /i = observed in |9j was caused by the 
overestimation of the nucleon mass Mm at large pion masses and the resulting overestimation of the coupling g. Now, 
with the more precise formula (13) for the nucleon mass, the transition is a crossover for all quark masses above the 
physical mass — 5 MeV. 



B. The deconfinement transition 



The deconfinement of the system is quantified by the Polyakov loop order parameter i and its conjugate f , which 
we treat as independent real variables. As with the chiral transition the transition temperatures could be determined 
locating the maxima of the temperature derivatives. However, as noted in our previous works [HI [9], there is in some 
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FIG. 2: The chiral phase diagrams of the models for several quark masses. The solid curves indicate a first order (discontinuous) 
transition with the critical endpoints marked by squares. Left: PLSM Right: PNJL 



cases a double peak structure in the derivatives with one peak coinciding with the chiral transition and caused by 
the interaction between the chiral and Polyakov sectors of the models. The second peak is a softer one and is related 
to the the transition present in the parametrization of the Polyakov loop mean field potential. Generally, the softer 
peaks occur at the values i — Xjl and I* — 1/2, so they could be considered as measures of the system becoming 
deconfined. This is a reasonable statement since confinement in the models is due to numerical suppression of quarks 
states through the order parameters I and i* . This means that alternatively the deconfinement could be considered to 
occur when the suppressing order parameters reach large enough value in order not to provide suppression anymore. 
Since the Polyakov loop order parameters £ and £* obtain, with our choice of potential, values roughly from to 1, 
the values £ = 1/2 and £* = 1/2 are a natural choice to be the indicator of when the system turns from a mostly 
confined state to a mostly deconfined state, bearing in mind that the transition is a crossover. To further simplify the 
analysis and readability of the figures we will, since £ and £* do not coincide at finite /i, define single deconfinement 
transition temperature as the average of the transition temperatures determined by £ = 1/2 and £* = 1/2. 

In Figures [3] and |4] the averaged deconfinement transition lines are shown for the physical quark mass « 5 McV 
and a larger quark mass niq = 50 MeV along with the corresponding chiral transitions. The first observation is that, 
defined this way, the deconfinement transition is independent of the chiral model and also of the amount of explicit 
chiral symmetry breaking i.e. the quark mass. This means that tuning the Polyakov potential, so that one obtains a 
coincidence of the deconfining and chiral transitions for one chiral model and a specific quark mass, will not give the 
same outcome in other cases. Furthermore, the deconfinement transition line, Tc^dccil^)j depends only very weakly 
on fi. This means that deconfinement and chiral restoration can be made to coincide only at /i = 0. However, on 
the basis of symmetries one does not expect deconfinement and chiral restoration to become independent at finite /i. 
Hence one is led to study possible /i-dependence in the Polyakov loop potential. 

To illustrate the effects of /x-dependence of the Polyakov potential on the deconfinement transition, we adopt 
the formulation from [T3] with the following modifications: In [T3] the /i-dependence of the Polyakov potential was 
through the critical temperature Tq, which could be described by the following parametrization 

Toifi) ^Tre-^^^°'°^^^''>\ (14) 

with the coefficient 6(/i) depending on the number of colors, massless flavors and the chemical potential. 

The parameters ao = 0.304 and Tr — 1.770 GeV were flxed to reproduce the Nf — lattice result TqIh = 0) = 270 
MeV. For two massless flavors this would mean To(/t = 0) = 208 MeV. However, since we want to make a comparison 
between the the /x-dependent and the /i~independent cases, we keep the /i = point as a reference point for the two 
cases. Hence we use 



(16) 
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FIG. 3: The phase diagrams of the models at the physical point m, ~ 5 MeV with no explicit /^-dependence in the Polyakov 
potential. The solid line is the chiral transition and the dashed line the deconfinement transition. The dotted lines correspond 
to quark number densities of 1 fm~^ and 2 fm"^. Left: PLSM Right: PNJL. 




FIG. 4: The phase diagrams of the models at niq = 50 MeV with no explicit /x-dependence in the Polyakov potential. The solid 
line is the chiral transition and the dashed line the deconfinement transition. The dotted lines correspond to quark number 
densities of 1 fm"^ and 2 frn'^ Left: PLSM Right: PNJL. 



This modification does not alter the fixing of the parameters ao and T^., which we fix to their above mentioned values. 
Our parametrization implies To(/.t = 0) = 270 MeV for any number of flavors. However, we are not interested in 
A^/-dependence since will exclusively consider the case Nf = 2 and we want to only use the jU-dependent Polyakov 
loop potential to illustrate the uncertainties which may arise in our consideration of (T, fj,)- phase diagrams. The 
A^/-dependence in (16) also neglects the effect of quark masses, which will suppresses Tq as discussed in [13]. This 
effect is not large for two light flavors, but grows more significant when considering larger quark masses. We stress 
that our intention here is simply to illustrate the effect of ^-dependent Polyakov potential on the phase diagram and 
for that purpose the simple parametrization we have chosen is sufficient. 

The phase diagrams obtained using the /x-dependent Polyakov potential are shown in Figures [5] and [6j The 
deconfinement transition now appears very different: The main new feature is that the deconfinement temperature 
follows the chiral restoration critical temperature more closely. Of course one may argue that this feature is put 
in by hand into the models of this type, but one the other hand the deconfinement phase transition is put in by 
hand already into the /x = Polyakov potential. Lattice determination of the coincidence of the chiral symmetry 
restoration and deconfinement at finite /i would provide strong motivation to use /^-dependent Polyakov potential in 
these effective modes. However, on the basis of the symmetries of the underlying gauge dynamics, one would indeed 
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expect deconfinement and chiral symmetry restoration to coincide also at finite /z. 



250 



200 



150 



100 



50 ■ 









s 

\. 

\ 

\. 

\ \ 
\ \ 
\ \ 
\ \ 
\ \ 


Const mass Tc 


\ \ 
\ 


Avg Polyakov Tc 


■ \ ^ 


Nq = 1 


■ \ \ 


Nq = 2/fm'^3 
1 1 1 1 ■ 1 1 1 1 





50 100 150 200 250 300 350 400 450 
\i [MeV] 



250 



200 



150 



100 



50 ■ 









\ \ 




\ \ 




\ \ 


Const mass Tc 


\ 


Avg Polyakov Tc 


\ 

\ ^ 


Nq = 1 


■ 1 ^ 


Nq = 2/fm'^3 
1 1 ■ 1 1 1 1 1 1 





50 100 150 200 250 300 350 400 450 
^l [MeV] 



FIG. 5: The phase diagrams of the models at the physical point rriq « 5 MeV with a /^-dependent Polyakov potential. The solid 
line is the chiral transition and the dashed line the deconfinement transition. The dotted lines correspond to quark number 
densities of 1 fm"^ and 2 fm"^. Left: PLSM Right: PNJL (Color online) 




FIG. 6: The phase diagrams of the models at = 50 MeV with a /i-dependent Polyakov potential. The solid line is the chiral 
transition and the dashed line the deconfinement transition. The dotted lines correspond to quark number densities of 1 fm~^ 
and 2 frn'^. Left: PLSM Right: PNJL (Color online) 

As a final application we will discuss tfie possibility of describing the quarkonic phase with these models. 



C. Quarkyonic matter? 

Recently there has been some interest towards a novel form of matter conjectured to exist in the QCD phase 
diagram jTU]. To briefly recall, the conjecture is based on considerations at 't Hooft large Nc limit and rests mainly 
on the following two features: First, since the free energy scales with the number of degrees of freedom, one has 
hierarchical contributions from the quarks and gluons, of the order of N,, and N^, respectively. Second, while the 
mesons and glueballs become free in the large Nc limit, their cubic and quartic interactions vanishing as 1/y/Wc and 
1/Nc respectively, the baryons remain strongly coupled. Both these features are, however, valid in the large Nc limit 
while for QCD Nc — 3. Furthermore, as we will now briefly discuss, the large Nc limit is not unique [21] . 

In the 't Hooft limit of large Nc the fermions are taken to transform according to the fundamental representation. 
This leads to the well known features at large Nc'. Planar diagrams dominate, and among these diagrams quark loops 
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are suppressed relative to gluonic ones. As the dimension of fermion representation is Nc while the one for gluons is 
— 1 ~ Nc, different orders of magnitude for the free energies emerge as Nc is taken large. 
However, consider the following redefinition of the quark fields 

Qa^^abcQbc^ Qbc^_Qcb^ 

i.e. consider quarks to transform in the two index anti-symmetric representation, which for SU(3) conceptually 
corresponds to renaming antiquarks as quarks. While this changes nothing in the dynamics for Nc = 3, the large 
Nc limit is entirely different. This is due to the fact that quarks in the anti-symmetric representation are counted 
similarly to gluons and hence they do not decouple at large Nc- The free energies of quarks and gluons are both of 
the order of Nc - The different behavior between mesons and baryons is similar to 't Hooft limit. 

These different large Nc limits emphasize different phenomenological features, and t is difficult to argue that either 
would be more realistic; more probably both are somewhat idealized and equally far from real three color QCD. And 
these two do not even exhaust the possible large Nc limits. For example, one can also consider the possibility of 
having a hybrid large Nc limit where out of three flavors one quark flavor transforms in the two-index antisymmetric 
representation while the other two are taken to transform in the fundamental representation. In this case it is possible 
to construct color singlet states consisting of three quarks at any Nc- These baryons are very different in comparison 
to baryons of the other two large Nc limits considered above: their masses do not grow with Nc and their Regge 
slopes coincide with the Regge slopes of mesons. 

We do not embark on a thorough analysis of the phenomenology associated with these large Nc limits; see e.g. j22j . 
From the above discussion we simply remark that the extrapolations from large Nc limits to real three color QCD 
should be taken with a grain of salt. A lattice study considering the quarkyonic phase in two-color QCD has recently 
appeared (23^. 

It has been suggested that that the quarkyonic phase can be characterized by a non-vanishing baryon density while 
the system is still in a confined phase |10j . In our models we have access to the quark number density 



an 



= 1^ (18) 



which we can use as a measure of the baryon density. 

If we first consider the case of a /i-independent Polyakov potential, illustrated in Figures [3] and |4j there is a sharp 
rise in the quark number density simultaneously with the chiral transition. As with the deconfinement transition, 
the absolute magnitude of this rise is not great, it is around 0.6 fm^"^, but the relative increase is about a factor of 
four. So it could be argued that at the point of the chiral transition also the number density changes from a nearly 
zero value to a non-zero one. Another way of characterizing this 'quarkyonic transition' is to assign a threshold value 
for Uq which separates the phases. This, however, is more arbitrary since the threshold value of riq one could choose 
is not in any way unique. In Figures [3] and [4] we have plotted the chiral transition line, which corresponds to the 
rise in quark number density, along with the deconfinement transition line, obtained as explained in the previous 
section, and two curves which correspond to values of 1 and 2 fm~^ of Hq. If one naively expects a baryon to have 
quark density around 3 fm~'^ then one should have nonzero net baryon density not later than when the average quark 
density hits 3 fm~'^, most likely even sooner. Seeing also that the deconfinement line lies at large temperatures, there 
is a substantial window for the quarkyonic matter to exist realizing the picture envisioned in [TU] . 

However, this picture changes considerably if the Polyakov loop potential depends explicitly on fi. This case is shown 
in Figures [5] and |6] First of all the quantitative behaviour of the quark number density changes quite drastically: Now 
Uq ~ 1 and Uq ~ 2 fm^"^ lines are well inside the chirally broken phase, but also the rise in quark number density 
associated with the chiral transition becomes significantly larger, in particular at large /i. But more importantly, 
the change in the deconfining transition, which now follows the chiral transition more closely, significantly decreases 
the area where possible quarkyonic matter could reside. Here one should remember that when implementing the 
/i-dependence to the Polyakov potential, we neglected effects from number of flavors and quark masses, which would 
bring the deconfinement lines down even faster than in Figures [5] and [6] and thus practically closing the window for 
the existence of quarkyonic matter. 



IV. CONCLUSIONS 



We have considered the (T, fi)- phase diagram of two-flavor QCD in effective models which take into account 
both chiral degrees of freedom relevant for the restoration of the chiral symmetry and Polyakov loop relevant for 
deconfienement. Earlier these studies have been performed by constraining these models to reproduce the physical 
vacuum, in the two flavor case essentially determined by the bare quark mass = m„ = (or alternatively by the 
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pion mass TOtt)- We have relaxed this assumption and treated niq as a free parameter of the model in order to study 
how the explicit chiral breaking manifests in the thermodynamics. 

We considered, side by side, two different models PNJL and PLSM which differ by the choice of the effective 
realization for the chiral sector. In earlier studies these two models have been shown to lead to qualitatively similar 
results for the thermodynamics both at finite temperature and density; the main quantitative difference has been 
shown to be in the location of the critical point in the (T, /i)- phase diagram. In this paper we have shown that, as 
a function of m^, a qualitative difference arises: while in PNJL model the critical point exists for any m^, in PLSM 
model the critical point exists only for a range of values of niq. In both models we find, for any m^, at most one critical 
point in contrast to the results in |20j where multiple critical points were observed for LSM model at non-physical 
values of the pion mass. 

Finally, we applied these effective models to consider the quarkyonic phase in two flavor QCD. The theoretical 
motivations on the existence of this novel phase are based on large A'c-limit of QCD and sensitive to which large Nc 
limit is considered model studies are required. Furthermore, the existence of this novel state of matter was shown to 
depend sensitively on the parametrization of the Polyakov loop potential at finite /i. In particular, for /i-independent 
Polyakov loop potential there appears a wide window in the (T, /i)- phase diagram for quarkyonic matter to exist 
while if fjL dependence in the Polyakov loop is introduced to obtain coincidence of chiral symmetry restoration and 
deconfinement similarly as at = 0, the window for quarkyonic matter practically closes. 
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